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Top quark mass measurement using the template method at CDF 
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We present a measurement of the top quark mass in the lepton+jets and dilepton channels of 
tt decays using the template method. The data sample corresponds to an integrated luminosity of 
5.6 fb _1 of pp collisions at Tevatron with y/s = 1.96 TeV, collected with the CDF II detector. The 
measurement is performed by constructing templates of three kinematic variables in the lepton+jets 
and two kinematic variables in the dilepton channel. The variables are two reconstructed top quark 
masses from different jets-to-quarks combinations and the invariant mass of two jets from the W 
decay in the lepton+jets channel, and a reconstructed top quark mass and m T 2 , a variable related to 
the transverse mass in events with two missing particles, in the dilepton channel. The simultaneous 
fit of the templates from signal and background events in the lepton+jets and dilepton channels to 
the data yields a measured top quark mass of M top = 172.1 ± 1.1 (stat) ± 0.9 (syst) GeV/c 2 . 

PACS numbers: 14. 65. Ha, 13.85.Ni, 13.85.Qk, 12.15.Ff 
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The top quark (t) is by far the heaviest known ele- 
mentary particle, with a mass approximately 40 times 
larger than the mass of its isospin partner, the bot- 
tom quark (6) [l]. The top quark contributes signifi- 
cantly to electroweak radiative corrections relating the 
top quark mass (M top ) and the W boson mass to the 
mass of the predicted Higgs boson within either the stan- 
dard model (SM) or beyond the SM di- Precision 
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measurements of M top provide therefore important con- 
straints on the Higgs boson mass in either model. Since 
the discovery of the top quark in 1995 |5| at the Fermilab 
Tevatron pp Collider, both the CDF and DO experiments 
have been improving the precision of the M top measure- 
ment [5] . However it is important to measure Mt op using 
different techniques and independent data samples in dif- 
ferent decay channels. Significant differences in the mea- 
surements of M top in different decay channels could indi- 
cate contributions from new physics beyond the SM [|]. 

This letter reports a measurement of the top quark 
mass using the template method We use samples 

of ti candidates in the lepton+jets and dilepton channels, 
corresponding to an integrated luminosity of 5.6 fb -1 of 
proton-antiproton collisions at -^/s = 1.96 TeV, collected 
by the CDF II detector 0- This is a general-purpose 
detector designed to study pp collisions at the Fermilab 
Tevatron. A charged-particle tracking system, consist- 
ing of a silicon microstrip tracker and a drift chamber, 
is immersed in a 1.4 T magnetic field. Electromagnetic 
and hadronic calorimeters surround the tracking system 
and measure particle energies. Drift chambers and scin- 
tillators, located outside the calorimeters, detect muon 
candidates. 



Assuming unitarity of the three-generation CKM ma- 
trix, the top quark decays almost exclusively into a W 
boson and a b quark [l[. The case where one W decays 
leptonically into an electron or a muon plus a neutrino 
and the other hadronically into a pair of jets defines the 
lepton+jets decay channel. The dilepton channel is de- 
fined as the case where both W’s decay leptonically into 
an electron or a muon plus a neutrino. 

Lepton+jets events are selected by requiring one iso- 
lated 0 electron (muon) with Et > 20 GeV (pr > 
20 GeV/c) and pseudorapidity |r;| < 1.1 [l2}. We also re- 
quire high missing transverse energy [H[ , > 20 GeV, 

and at least four jets. Jets are rec onstructed with a cone 
algorithm 14 with radius R = \J (A??) 2 + (A0) 2 = 0.4. 
Jets originating from b quarks are identified (tagged) 
using a secondary vertex tagging algorithm [15(. We 
request at least one jet to be tagged as a b jet. We 
divide the sample of candidate lepton+jets events into 
sub-samples of one 6-tagged jet (1-tag) and two or more 
6-tagged jets (2-tag). In events with more than two b- 
tagged jets, we consider the two highest Et jets as b 
quark candidates and treat the other 6-tagged jets as 
non 5-tagged jets. In the 1-tag sample, we require ex- 
actly four jets with transverse energy Et > 20 GeV and 
\r/\ < 2.0. In the 2-tag sample, three jets are required to 
have Et > 20 GeV and \r/\ < 2.0, and at least one more 
jet is required to have Et > 12 GeV and |?y| < 2.4. We 
apply an additional cut on the scalar sum of transverse 



energies in the event, Ht = E l ^ pton 



$T + J2 



E 



+t 



jets ^ T ’ 



requiring Ht > 250 GeV for all events to further reject 
backgrounds. E™ uon = p™ uon is assumed in the Ht cal- 
culation. 



The primary sources of background in the lepton+jets 
channel are W+jets and QCD multijet production. We 



also consider small contributions from Z-\- jets, diboson, 
and single-top production. To estimate the contribution 
of each process, we use a combination of data and Monte 
Carlo (MC) based techniques described in Ref. [16]. For 
the Z-\- jets, diboson, single top, and tt events we nor- 
malize MC simulation events using their respective theo- 
retical cross sections 0®. QCD multijet background 
is estimated using the data referring to techniques de- 
scribed in Ref. 20]. The shape of W+jets background is 
obtained from MC while the number of W+jets events 
is determined from the data by subtracting all the other 
backgrounds and tt. 



Three observables are used from each lepton+jets 
event: two reconstructed top quark masses ( m\ eco and 
TO( CC °( 2 )) and the invariant mass of the two jets from the 
hadronically decaying W boson ( rrijj ). We have complete 
reconstruction of the ti kinematics in the lepton+jets 
channel m with constraints from the precise W bo- 
son mass and requiring the t and i masses to be the 
same. With the assumption that the leading four jets 
in the detector come from the ti decay products, there 
are six and two possible assignments of jets to quarks 
for 1-tag and 2-tag respectively. A minimization is per- 
formed for each assignment using a y 2 comparison to the 
ti hypothesis with m\ eco taken from the assignment that 
yields the lowest y 2 . To increase the statistical power 
of the measurement, we employ an additional observ- 
able rri t eco< ' 2 ^ from the assignment that yields the 2 nd 
lowest x 2 - Events with the lowest y 2 > 9.0 are re- 
moved from the sample to reject poorly reconstructed 
events. The dijet mass rrijj is calculated as the invari- 
ant mass of two non 6-tagged jets which provides the 
closest value to the world average W boson mass of 
80.40 GeV/c 2 [J. We apply boundary cuts on m r t eco and 
m T t eco(2) (100 GeV/c 2 < m\ eco ,rn t eco(2) < 350 GeV/c 2 ) 
and rrijj (50 GeV/c 2 < rrijj < 120 GeV/c 2 for 1-tag 
events and 50 GeV/c 2 < rrijj < 125 GeV/c 2 for 2-tag 
events), and normalize the probability density function in 
the signal region. The estimated number of background 
events and the observed number of events after event se- 
lection, x 2 cu t) and boundary cuts are listed in Table [Q 
for the lepton+jets decay channel. 



To select dilepton candidate events, we require two op- 
positely charged leptons with Et > 20 GeV (for elec- 
trons) or pt > 20 GeV/c (for muons). One lepton is 
required to be isolated in the central region (|? 7 | < 1.1) of 
the detector, but the other can be a non-isolated lepton 
in the central region or an isolated electron in the forward 
region (1.1 < |? 7 | < 2.0). We also require #t > 25 GeV, 
and at least two jets with Et > 15 GeV and |?y| < 2.5. To 
further reject backgrounds, we require Ht > 200 GeV. 
In measuring the top quark mass, we divide the dilep- 
ton sample into events with 6-tagged jets (tagged) and 
without 6-taggedjets (non-tagged) . 

Drell-Yan, diboson, and W+jets (fake lepton) events 
are the primary sources of background in the dilepton 
channel. We estimate the rate of the Drell-Yan and di- 
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TABLE I: Expected and observed numbers of signal and 

background events assuming tt production cross section a t j = 
7.4 pb and Mt op = 172.5 GeV/c 2 in the lepton+jets channel. 





1-tag 


2-tag 


JF+jets 


53.4 ± 17.5 


8.5 ± 3.0 


QCD multijet 


13.1 ± 10.6 


1.8 ± 1.5 


Z + jets 


4.7 ± 1.0 


0.5 ± 0.1 


Diboson 


6.3 ± 0.8 


0.8 ± 0.1 


Single top 


4.9 ± 0.4 


2.0 ± 0.2 


Background 


105 ± 21 


14.2 ± 3.3 


tt signal 


590 ± 74 


293 ± 45 


Expected 


694 ± 77 


307 ± 45 


Observed 


695 


286 



TABLE II: Expected and observed number of signal and 

background events assuming tt production cross section ojj 
= 7.4 pb and M top = 172.5 GeV/c 2 in the dilepton channel. 




non-tagged 


tagged 


Diboson 

Drell-Yan 

LF+jets (fake lepton) 


19.2 ± 3.3 
31.5 ± 3.9 
30.8 ± 9.4 


0.7 ± 0.2 
3.7 ± 0.2 
4.6 ± 1.3 


Background 
tt signal 


81.6 ± 10.4 
124 ± 16 


8.9 ± 1.4 
151 ± 19 


Expected 

Observed 


205 ± 19 
237 


160 ± 19 
155 



boson events with calculations based on MC simulations. 
For the Droll- Yan Z+jets process, we normalize the MC 
sample by matching the number of Z events predicted 
and observed in the Z mass region between 76 GeV/c 2 
and 106 GeV/c 2 . We use data to estimate the rate of 
W+jets (fake lepton) events where an event has one real 
lepton and one of the jets misidentified as the other lep- 
ton. The detailed procedure of background estimation in 
the dilepton channel is described in Ref. [2lJ. For each 
event we calculate a reconstructed top quark mass m^ WA 
using the neutrino weighting algorithm [22} , and we cal- 
culate a quantity tot 2 [23] . Here m/p 2 is a variable re- 
lated to the transverse mass of the mother particles (top 
quark in the tt system) in events with two missing par- 
ticles from pair production of the mother particles. We 
firstly use this variable for the top quark mass measure- 
ment in the dilepton channel 9]. We require these ob- 
servables to be consistent with the top quark signal by 
demanding 100 GeV/c 2 < mf WA < 350 GeV/c 2 and 
30 GeV/c 2 < TOp 2 < 200 GeV/c 2 . The estimated num- 
ber of background events and the observed number of 
events after event selection are listed in Table m for the 
dilepton decay channel. 

We estimate the probability density functions (p.d.f.’s) 
of signal and background using kernel density estima- 
tion (KDE) jH] that constructs the p.d.f. without any 



assumption of a functional form. In the lepton+jets chan- 
nel, we use the three dimensional KDE that accounts for 
the correlation between the three observables. In the 
dilepton channel, instead, we use the two dimensional 
KDE. The dijet mass rrijj of the two jets assigned to 
the W in the lepton+jets channel is used for in situ cal- 
ibration of jet energy scale (JES) 0, Q. The p.d.f.’s 
for the observables are estimated at discrete values of 
M t0 p from 130 GeV/c 2 to 220 GeV/c 2 , with increments 
from 0.5 GeV/c 2 in the region immediately above and 
below 172.5 GeV/c 2 to 5 GeV/c 2 near the extreme mass 
values, and at discrete values of Ajes from —3.0 a c to 
3.0 <t c with increments of 0.2 0 + where cr c is the CDF 
JES fractional uncertainty [25} and Ajes corresponds to 
the difference between the energy scale in MC simula- 
tion and data. We interpolate the MC distributions to 
find p.d.f.’s for arbitrary values of M top and Ajes us- 
ing the local polynomial smoothing method [26}. We fit 
the signal and background p.d.f.’s to the distributions of 
the observables in the data using an unbinned maximum 
likelihood fit [27} where we minimize the negative loga- 
rithm of the likelihood using MINUIT [28| . The likelihood 
is built for each sub-sample separately, 1-tag and 2-tag 
for lepton+jets events, non-tagged and tagged for dilep- 
ton events, and an overall likelihood is then obtained by 
multiplying them together. We independently obtain the 
results from the lepton+jets channel, the dilepton chan- 
nel, and the two channels combined. In the combined fit, 
the dilepton channel uses the JES calibration found in 
the lepton+jets channel. We evaluate the statistical un- 
certainty on M top by searching for the points where the 
negative logarithm of the likelihood exeeds the minimum 
by 0.5. Ref. 0 , 0 provides detailed information about 
this technique. 

We test the mass fit procedures using 3000 pseudoex- 
periments for a set of 14 different M top values ranging 
from 159 GeV/c 2 to 185 GeV/c 2 . In each experiment, 
we select the number of background events from a Pois- 
son distribution with a mean equal to the expected total 
number of background events in the sample and the num- 
ber of signal events from a Poisson distribution with a 
mean equal to the expected number of signal events nor- 
malized to a tt pair production cross section of 7.4 pb at 
M top = 172.5 GeV/c 2 [h|. The distributions of the av- 
erage mass residual (deviation from the input top mass) 
and the width of the pull (the ratio of the residual to the 
uncertainty reported by MINUIt) for simulated experi- 
ments are corrected to be unity and zero respectively. 
The corrections are m corr = 1.04 x m m eas ~ 6.8 GeV/c 2 , 
m cor r = 1.03 x mmeas — 5.5 GeV/c 2 , and m corr = 
1.03 x mmeas ~ 5.9 G eV/c 2 for combined fit, lepton+jets, 
and dilepton channel respectively, where m meas is the 
raw value from likelihood fit and m corr is the corrected 
value of the measurement. We increase the measured un- 
certainty by 4% for combined fit and lepton+jets channel 
and 3% for dilepton channel to correct the width of the 
pull. 

We examine various sources of systematic uncertainties 



6 







FIG. 1: Distributions of the three variables used to measure Mt op in the lepton+jets channel, showing 1-tag and 2-tag samples 
separately. The data are overlaid with the predictions from the KDE probability distributions using Mt op = 172.0 GeV/c 2 and 
the full background model. 



that could affect the measurement by comparing the re- 
sults of pseudoexperiments in which we vary relevant pa- 
rameters within their systematic uncertainties. The dom- 
inant sources of systematic uncertainty are the residual 
JES d,[25j] and signal modeling. We vary JES parameters 
within their uncertainties in both signal and background 
MC generated events and interpret the shifts in the re- 
sults of the pseudoexperiments as uncertainties. For the 
dilepton channel, which has no in situ calibration, the 
JES is the single dominant uncertainty. The uncertainty 
arising from the choice of MC generator (signal model- 



ing) is estimated by comparing the results of pseudoex- 
periments generated with PYTHIA [29| and HERWIG [30]. 
The b - JES systematic uncertainty arising from our mod- 
eling of b fragmentation, b hadron branching fractions, 
and calorimeter response captures the additional uncer- 
tainty not taken into account in the (residual) JES. We 
estimate the systematic uncertainty due to modeling of 
initial-state gluon radiation and final-state gluon radia- 
tion by extrapolating uncertainties in the pr of Drell- 
Yan events to the tt mass region Q. We estimate the 
systematic uncertainty due to parton distribution func- 
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FIG. 2: Distributions of the two variables used to measure M top in the dilepton channel, showing non-tagged and tagged samples 
separately. The data are overlaid with the predictions from the KDE probability distributions using Mt op = 170.0 GeV/c 2 and 
the full background model. 



TABLE III: Estimated systematic uncertainties in the com- 
bined fit (Comb), lepton+jets (LJ), and dilepton (DIL) (unit 
in GeV/c 2 ). 



Source 


Comb 


LJ 


DIL 


(Residual) Jet Energy Scale 


0.5 


0.5 


3.0 


Signal modeling 


0.7 


0.7 


0.3 


b Jet energy scale 


0.3 


0.3 


0.4 


Initial and final state radiation 


0.1 


0.1 


0.2 


Parton distribution functions 


0.1 


0.1 


0.3 


Gluon fusion fraction 


<0.1 


<0.1 


0.1 


Lepton energy 


<0.1 


0.1 


0.3 


Background shape 


0.1 


0.1 


0.3 


Multiple hadron interaction 


0.2 


0.1 


0.2 


Color reconnection 


0.2 


0.2 


0.5 


MC statistics 


0.1 


0.1 


0.1 


Total systematic uncertainty 


0.9 


0.9 


3.1 



tions (PDFs) by varying the independent eigenvectors of 
the CTEQ6M [3l| PDFs, varying A pen, and compa ring 
the results obtained with CTEQ5L [32j and mrst72 |33| 
PDFs. In estimating the systematic uncertainty from the 
top quark production mechanism, we vary the fraction 



of top quarks produced by gluon-gluon annihilation from 
6% to 20%, corresponding to the one standard deviation 
upper bound on the gluon fusion fraction [34j. We esti- 
mate systematic uncertainties due to the lepton energy 
and momentum scales by propagating shifts in electron 
energy and muon momentum scales within their uncer- 
tainties. Background shape systematic uncertainties ac- 
count for the variation of the background composition. 
We estimate the multiple hadron interaction systematic 
uncertainty to account the effect from the difference in 
the average number of interactions between MC samples 
and the data. The color reconnection (CR) systematic 
uncertainty [§j| is evaluated by MC samples generated 
with and without CR effects using different tunes [36] of 
PYTHIA. The total systematic uncertainties, adding indi- 
vidual components in quadrature, are 0.9 GeV/c 2 in the 
combined fit, 0.9 GeV/c 2 in the lepton+jets channel, and 
3.1 GeV/c 2 in the dilepton channel. 

We perform the likelihood fits to the data using the 
observables discussed in this letter and apply the cor- 
rections obtained using the simulated experiments. We 
obtain for the lepton+jets channel, a top quark mass 

M top = 172.2 ± 1.2 (stat) ± 0.9 (syst) GeV/c 2 
= 172.2 + 1.5 GeV/c 2 , 



while for the dilepton channel, 

M top = 170.3 ± 2.0 (stat) ± 3.1 (syst) GeV/c 2 
= 170.3 ±3.7 GeV/c 2 . 

The two channel combined fit yields a top quark mass 

M top = 172.1 ± 1.1 (stat) ± 0.9 (syst) GeV/c 2 
= 172.1 ± 1.4 GeV/c 2 . 

Figure Q] shows the measured distributions of the observ- 
ables used for the M top measurement in the lepton±jets 
channel overlaid with density estimates using ti signal 
events with M top = 172 GeV/c 2 (close to the measured 
M top hr the lepton±jets channel) and the full background 
model. Figure [2] shows the corresponding distributions in 
the dilepton channel using ti signal events with M top = 
170 GeV /c 2 (close to the measured M top in the dilepton 
channel) . 

In conclusion, we have performed a measurement of 
the top quark mass using the template method simulta- 
neously in the lepton±jets and dilepton channels. The 
result, M top = 172.1 ± 1.4 GeV/c 2 , is consistent with 
the most recent world average of M top = 173.3 ± 
1.1 GeV/c 2 [5]]. In the lepton+jets channel, we use the 
same data set as the best single M top measurement [37], 
and have a consistent result with slightly larger uncer- 
tainty. In the dilepton channel, we achieve the single 



most precise measurement of M t0 p in this channel to date 
and the result is in good agreement with the measure- 
ment in the lepton±jets channel. 



Acknowledgments 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu- 
tions. This work was supported by the U.S. Depart- 
ment of Energy and National Science Foundation; the 
Italian Istituto Nazionale di Fisica Nucleare; the Min- 
istry of Education, Culture, Sports, Science and Tech- 
nology of Japan; the Natural Sciences and Engineering 
Research Council of Canada; the National Science Coun- 
cil of the Republic of China; the Swiss National Science 
Foundation; the A.P. Sloan Foundation; the Bundesmin- 
isterium fur Bildung und Forsclrung, Germany; the Ko- 
rean World Class University Program, the National Re- 
search Foundation of Korea; the Science and Technology 
Facilities Council and the Royal Society, UK; the Institut 
National de Physique Nucleaire et Physique des Partic- 
ules/CNRS; the Russian Foundation for Basic Research; 
the Ministerio de Ciencia e Innovation, and Programa 
Consolider-Ingenio 2010, Spain; the Slovak R&D Agency; 
the Academy of Finland; and the Australian Research 
Council (ARC). 



Ill C. Amsler et al. (Particle Data Group), Phvs. Lett. B 
667 , 1 (2008). 

[2] ALEPH, CDF, DO, DELPHI, L3, OPAL, SLD, the LEP 
Electroweak Working Group, the Tevatron Electroweak 
Working Group, and the SLD Electroweak and Heavy 
Flavor Working Groups, arXiv:1012.2367v2. 

[3] H. Flacher et al, Eur. Phys. J. C 60 , 543 (2009). 

[4] F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 
74 , 2626 (1995); S. Abachi et al. (DO Collaboration), 
Phys. Rev. Lett. 74 , 2632 (1995). 

[5] The Tevatron Electroweak Working Group (CDF 
and DO Collaborations), FERMILAB-TM-2466-E, 
arXiv:1007.3178vl. 

[6] G. L. Kane and S. Mrenna, Phys. Rev. Lett. 77 , 3502 
(1996). 

[7] A. Abulencia et al. (CDF Collaboration), Phys. Rev. D 
73 , 032003 (2006). 

[8] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. D 
79 , 092005 (2009). 

[9] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. D 
81 , 031102 (2010). 

[10] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71 , 
032001 (2005). 

[11] A lepton is isolated when [p!]? tal — p!^ pton J /p^ pton < 

0.1, where py tal is the total transverse momentum (en- 
ergy) and py pton is the lepton transverse momentum (en- 
ergy) for muon (electron) in a cone of radius A R = 
( A77) 2 + (A^) 2 = 0.4 [12] with axis along the direc- 



tion of the lepton. 

[12] We use a right-handed spherical coordinate system with 
the origin at the center of the detector with the z-axis 
along the proton beam and the y-axis poiting up. 9 
and 4> are the polar and azimuthal angles, respectively. 
The pseudorapidity is defined by 77 = — lntan(#/2). The 
transverse momentum and energy of a detected particle 
or jet are defined by pr = psin($) and Et = E sin(0), re- 
spectively, where p and E are the momentum and energy 
of the particle. For the reconstructed top quark decay 
products used in the mx 2 calculation, the transverse en- 
ergy is defined by Et = \/rn 2 + p y, where m is the mass 
of the product. 

[13] The missing transverse energy, an imbalance of energy in 
the transverse plane of the detector, is defined by $t = 

Stowers |, where hr is the unit vector normal to 

the beam and pointing to a given calorimeter tower and 
Et is the transverse energy measured in that tower. 

[14] F. Abe et al. (CDF Collaboration), Phys. Rev. D 45, 
1448 (1992). 

[15] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71 , 
052003 (2005). 

[16] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 
105, 012001 (2010). 

[17] J. M. Campbell and R. K. Ellis, Phys. Rev. D 60, 113006 
(1999). 

[18] B. W. Harris el al., Phys. Rev. D 66, 054024 (2002). 

[19] S. Moch and P. Uwer, Nucl. Phys. B, Proc. Suppl. 183, 
75 (2008). 



9 



[20] T. Aaltonen el al. (CDF Collaboration), Phys. Rev. D 
77 , 011108 (2008). 

[21] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. D 
82 , 052002 (2010). 

[22] B. Abbott et al. (DO Collaboration), Pliys. Rev. D 60, 
052001 (1999); A. Abulencia et al., Phys. Rev. D 73, 
112006 (2006). 

[23] C. Lester and D. Summers, Phys. Lett. B 463, 99 (1999); 
A. Barr, C. Lester, and P. Stephens, J. Phys. G 29 , 2343 
(2003). 

[24] K. Cranmer, Comput. Phys. Commun. 136 198 (2001). 

[25] A. Bhatti et al, Nucl. Instrum. Methods Phys. Res. Sect. 
A 566 375 (2006). 

[261 C. Loader, Local Reqression and Likelihood (Springer, 
New York, 1999). 

[27] R. Barlow, Nucl. Instrum. Methods Phys. Res. Sect. A 
297 496 (1990). 



[28] F. James and M. Roos, Comput. Phys. Commun. 10 343 
(1975). 

[291 T. Sjostrand, S. Mrenna, and P. Skands, J. High Energy 
Phys. 05 (2006) 026. 

[30] G. Corcella et al., J. High Energy Phys. 01 (2001) 010. 

[31] J. Pumplin et al., J. High Energy Phys. 07 (2002) 012. 

[32] H. L. Lai et al. (CTEQ Collaboration), Eur. Phys. J. C 
12 , 375 (2000). 

[33] A. D. Martin et al., Eur. Phys. J. C 14 , 133 (2000). 

[34] M. Cacciari et al., J. High Energy Phys. 04 (2004) 068. 

[35] P. Z. Skands and D. Wicke, Eur. Phys. J. C 52, 133 
(2007). 

[36] P. Z. Skands, Phys. Rev. D 82 , 074018 (2010). 

[37] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 
105 , 252001 (2010). 




